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ABSTRACT: 

Theoretical expressions are derived for the sound scattered 
at a time-dependent rough surface. The calculations are made for 
a Gaussian shaded source transducer and point receiver. The 
Helmholtz theorem and Fresnel approximation are used. The rough 
surface is assumed to be a traveling wave and to have a traveling 
wave packet type of correlation function. The coherent component 
of the signal is the product of the Fourier transformation of the 
surface distribution function and the smooth surface reflection 
signal. Comparison of theory and experiment shows the coherent 
component to be sensitive to the non-Gaussian character of the 
wind-blown water waves. The incoherent components and the temporal 
correlation function of the scattered sound are given. For the 
special case of a traveling cosine wave type of rough surface, the 
spectrum of the scattered sound includes components which are mul- 
tiples of the frequency of the surface wave. For surfaces describ- 
able by a bivariate Gaussian distribution function, the temporal 
correlation is a function of, but not the same as, the time de- 
pendence of the rough surface. The scattered sound is insensitive 
to the spatial correlation function of the surface at distances 
larger than the dimensions of the transducer divided by the cosine 
of the incident angle. The final expressions are complex error 
integrals and can be used for all values of roughness. This task 
was supported by Naval Ship Systems Command (Code PMS 388) . 
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I • EXPERIMENT AND THE INVERSE PROBLEM 



Sound scattered at a time dependent surface takes on some of the 
time dependence of the surface. For example, the upward or downward motion 
of the surface causes the frequency of the reflected signal to be Doppler 
shifted. It is easy to demonstrate the effects in laboratory experiments 
(Fig. 1.1). Small waves on the water cause the phases of the signals to 
fluctuate. Experiments have shown the reflected signals are (crudely 
stated) modulated by the surface and often the spectrum 
of the surface can be identified as a component of the spectrum of the 
reflected signal. With our growing technology in remote sensing, the 
importance of being able to go from scattering measurements back to a 
description of the time dependent surface is extremely important. It 
may be more important to know when and how far one can go for a given 
technique . 

Doing the inverse problem requires a complete knowledge of the relation 

ship of the measurement to the quantity being sensed. Too often, the 

result of having this knowledge is "I can't measure what I want by doing 

the experiment I am doing." For example? in the plane wave approximation, 

Eckart (1953) showed that mean square scattered signals are simply related 

to the correlation function of the surface at small ycr 

(y = k(cos Q^+ cos9 2 )/2, k is 2^/wavelength, o is rms roughness, and 

9^ and 9^ are the incident and reflected angles). Unfortunately the 

2 2 

scattered signal is proportional to y a and correspondingly very small. 

If the source and receiver are omnidirectional or broad beamed, the 
scattered sound is usually identified as being a fluctuation of the signal 
and is often inseparable from the noise. 



4 



In this paper I add another limitation (which is derived in the 
Fresnel approximation. The scattered sound is insensitive to spatial 
correlation function of the rough surface at correlation distances larger 
than the dimensions of the transducer/cos9^ . 

I believe a bit of discussion of the theoretical problem and 
approximations is needed. In seeking the source of a discrepancy between 
my theory (1971) and some of our laboratory measurements, G. A. Sandness 
identified the calculation of the incident sound at the surface as being 
the difficulty for two reasons, (private communication). First, in using 
the Helmholtz equation, the incident sound signal should satisfy the wave 
equation. The combination of a point source and directional function (or 
illumination function on the surface) does not satisfy the wave equation. 
Second, the limitation of the expansion to second order terms may be a 
very poor approximation for a diverging wave at the surface. 

Another way of thinking about this is to regard the phase and amplitude 
of the incident signal as being the hologram of the source. The hologram 
of a finite object is different from that of a point source. The reader 
may wish to read Melton and Horton (1970). 

Having said that there are difficulties with the Fresnel approximation, 
I have chosen to use it. Also, I use a Gaussian shaded source transducer 
and a point receiver. This shading fits the measured response of our 
transducers and facilitates repeated integrations. I have obtained ex- 
pressions that can be used over a wide range of surface roughnesses and 
correlation functions without resorting to separate expansions for the high 
and low frequency limits. The end result is the mean square scattered sound 
and its temporal correlation function. 
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II. SCATTERING FUNCTION FROM THE HELMHOLTZ EQUATION 



The development of a theoretical relationship of a scattering function 
to the properties of the surface requires analysis of the scattering of 
acoustic signals by rough surfaces. This has been the subject of a number 
of papers and is treated in several books, notably Beckmann and Spizzichino 
(1963), Tolstoy and Clay (1966), Ol'Shevskii (1967), Fortuin (1970), and 
Horton (1971). Although we will not discuss it here, in the electronic 
and radio engineering journals an extensive literature exists on the 
scattering of electromagnetic waves by various types of irregular and rough 
surfaces. We are obliged to present the theory in some detail because the 
relevant underwater acoustics theory has not been applied to these types 
of problems. Most of the studies have been empirical. 

The derivation of the scattering equation from the Helmholtz theorem 
is usually based upon the assumption of local plane reflections. It is 
suggested that the reader who is interested in the details refer to the 
development given by Tolstoy and Clay (1966), who base their derivations 
on those of Eckart (1953) and Beckmann and Spizzichino (1963). The sound 
absorption will be ignored in the derivation. 

The general assumptions are: 1) The source and receiver are far from 

the illuminated area. 2) The dimensions of the source are small compared 
to R^ and R^. 3) The source is a Gaussian shaded transducer and is directed 

along R^, (Fig. 2.1). Also, the individual elements are in phase and 
incident sound pressure at (x,y,£) is the integral over ds 1 . 4) No shadows 

are present. 
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Fig. 2.1 Geometry 

ds is at the position x, y, £ relative to the origin. The 
plane of the source is perpendicular to R^. 
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Because the surface is rough, we can expect that at grazing angles 
some of the surface will be in shadow. The proportion of the area in 
shadow to the total area scanned depends upon the shape of the surface 
and on the grazing angle. Wagner (1967) discussed this for rays and a 
randomly rough surfaces. In any case, the assumption of no shadowing 
is violated. It is also likely that the radii of curvature of some 
features on the surface are small compared to the acoustical wave- 
length. In that case, the sound does not have a local plane reflection 
and in order to describe the scattered sound we must use higher order 
expansions. 

Eckart (1953) remarked that the boundary conditions are troublesome. 
For example, one can set' 



P = ^P 1 




( 2 . 1 ) 

( 2 . 2 ) 



where p^ and p are the incident and reflected signals on the local surface 
and d/ is the normal derivative. Using ft = -1 for a free surface, Eckart 
suggested that in the smoother areas, the second condition might hold. In 
deep shadows, he thought that 

= (2 - 3) 



j2& 

Sn 



might be reasonable. Horton and Muir (1967) assumed an average boundary 
condition and combined the two equations involving the normal derivatives 



to obtain 



ft = 



Po Co cos 0- - P-C.COS0 

2 2 1 11 r 

P 0 C COS 0 . + P-C-COS 0 

2 2 1 11 r 



reflection coefficient for plane waves. 



here 1 and 2 refer to medium 1 and 2. 
in medium 2. 



0 is refracted direction 
r 
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(2.4) 



c^/dn = 0 or <$p/bn = 0 

They also used Eckart’s small slope approximation and replaced the 
normal derivative by 5/3z. Horton et al . (1967) compared theoretical 
computations and the experimental scattered sound and found the 
agreement to be quite good. The slope correction can be included 
by performing an integration by parts [Tolstoy and Clay, p. 196-199 
(1966)]. In the specular direction, the result is the same as 
obtained with S/Sn ~ S/Sz. The result is different for backscattered 
sound . 

The most direct way to demonstrate the approximations is to 
start with the integral expression for the scattered acoustical 
pressure. With the aid of the Helmholtz theorem [Born and Wolf, 
p. 375 (1965)] it is 



where 




B 2 = JIpc (2n) _1 
II = source power 



(2.5) 

( 2 . 6 ) 

(2.7) 



and the normal is drawn toward the receiver. R and R are the source 
and receiver distances to ds (at x, y, z) . In applying the Helmholtz 
theorem, the surface is closed by a hemisphere at infinity. There 
are no other sources and all waves are outgoing from our source. 

Hence the contribution to the integral for the surface at infinity can 
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be ignored and the integral over the illuminated interface is used. 
We choose to use the less conventional form of the normal toward 
the receiver because that is the direction of the wave propagation 
and it does not matter where the source is placed. The substitution 
of (2.1) and (2.2) into (2.5) gives 



■ zk-f* lr (p i u) l s ds 



( 2 . 8 ) 



or the Horton-Muir condition, 2.4 gives 

1 P ri c>U I 

p ■ urj s * p i xi s ds 



(2.9) 



After expansion of the integrands for a moderately directive source, 
the integrals will be the same and the kind of boundary condition 
depends upon a constant factor in front of the integral. 

There are several ways of writing the algebra for the expansion 



and they are all messy. Since the transducer plane is perpendicular 

r 



' 2 2 

to R , R and R are 



*2 1 2 1 2 
R = (-R^ s ^ n ®i + x cos ~ + (y - y) + (R^cos0^ - 

R 2 = (R sin© cos© - x; 2 -b (R sin© sin© - y) 2 + (R cos© - (; 
223 223 2 1 



( 2 . 10 ) 



With the aid of the binomial expansion and retaining £ and second 

i 

order terms, R and R are 

. '2, '2 x 2 cos 2 © 1 + y 2 xx cos©, ' 

n ' D , * + V I L , . l . TJL 

- 1 2R l 2R R x R l 



+ x sin©^ - £ cos©^ + . . . 



x 



R ~ R„ + (l - sin 2 © 2 cos 2 © 3 ') + (l - 



sin © 2 sin © 



2 ~ ) 
3' 



- x sin© 2 cos© 3 - y sin© 2 sin© 3 - £ cos© 2 + . . . 
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x sin©^-£ 
2 



( 2 . 11 ) 



For small slopes, the normal derivative is approximately (or 5/ 

here). On making this approximation and also including the integral over 



the transducer ds 1 , one obtains 



p:- 



ik B$F e 



-iO)t 



2K R 1 R 2 



P , , f ik (R ' + R) 

> ds j e 



ds 



(2.12) 



(cos0^ + cosG^)/2, boundary cond. (2.1 and 2), d/ dn ~ d/dz 



_ / 1 + co s 0- + cos0 o - sin0- sin0 o cos0 o 1 . . /0 . oX 

F = < 1 2 1 2 3 , slope correction (2.13) 



cos£^ + cosO^ 



(Tolstoy and Clay, 1966) 



cosG^/ly boundary cond. (2.4) 



g = ^ exp 



x 



2 2 
i ^ __ i 



(2.14) 



W 



Integration over the source yields (after algebraic manipulation) 

fiat -ikC^ + R 2 )J 



p ~ - 



ik B£RF e 



2” W 1 - i d L > ! a - 



exp a^x 2 - a^y 2 + 2 iccx + 2 ipy + 2 iy^j dydx 

ik cos 2 0 1 R fl - i (1 - R„/R ) d T 1 
_ 1 x L 2 x L J 



a x - ' 2 R 1 R 2 



1 - id T 



ik R v [i - i (i - R 2 /yd w ] 



a y ~ " 2 R, R 
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1 - id, 



W 



(2.15) 



= R 2 + R-^ (sin 2 © 2 cos 2 0g)/cos 2 0^ 
Ry = R 2 + R^ (1 - sin 2 0 2 sin 2 ^) 
d L = kL 2 /(2R 1 ) d w s kW 2 /(2R 1 ) 



( 2 . 16 ) 



2a = k(sin6^ - sin^ cos ®j) 
2(3 = -k(sin02 sin0^) 



2y = -k(cos9^ + cos^) 
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Integration of 2.15 for £ = 0 yields the image solution. Although 
the size of the transducer is included and the expressions are dependent 
on L and W, I don ! t think the expressions are accurate for near 
field computations because the higher powers of x and y were dropped 
in 2.11. 

Proceeding, the covariance of the sound pressure is 



(t + T) p" (t)^ = 



. 2 2 2 2 -io)r 
= k B F e 



2 2 2 2 3 * 2 1 ' 

4 * \\ (1 + d L Z )*Q + d w V 



2 * , 2 2 * i2 

-ax-ax'^-ay - ay 1 

xx y y 



+ 2ia(x-x') + 2iP(y-y') + 2i.y(£ -£') J 

dy dx dy' dx'^ (2.17) 

I assume the surface is random and £ has a zero mean. The correlation 



function of £ is 



X Y T 
2 2 2 



a \ |r = lim 

X,Y,T 



sin 



£ dx dy dt 



(2.18) 



X Y T 
"2 2 2 



where in functional notation ^ could be written as t,). 

Assuming a Gaussian surface, the bivariate Gaussian probability density is 

2-1 2 - 3 - 

W = (2na ) (1-t ) 2 

exp | - [ 2 (l-* 2 )a 2 ] _1 [ £ 2 + £' 2 - 2^'+] 3 (2.19) 



Since the only random quantities are £ and the average in (2.17) 

operates on them as follows: 



For a given roughness wave length A, the shorter wave length roughness 
in an area having the dimensions of several A has the same statistical pro- 
perties as all other similar size areas. 
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/ e 2i r<?-5'>) ■ fj He 2ir(E-C) 



d£d£' °2 



= exp [ - 4r 2 0 2 (l-'IO] 



< > 

The following changes of variable permit integration over x" and y". 



( 2 . 20 ) 



x = x" + 1/2 



x’ = x” - i/2 



y - y" + t}/2 



( 2.21 ) 



y '=* y" - n /2 

Jio substitution of 2.20 and 2.21 into 2.17 and evaluation of the 



integral yields, after the usual algebra, 
^p(t + t) p (t)^ ~ 



where 



b 2 fV. ^ 

2jtR^ LW cos0^ 



/ [• v 2 -v 2 + 2ia|+ 2ipT! - 4r 2 a 2 (l-il/)j dl dl} (2.22) 



a fc = 



2 2 
R cos0 
x 



1 E * d^d-R;/!^) 2 ] 



5 .. 2 2 

2 R 2 L 

R 2 r 1 + d 2 (1-R 2 /r v ) 2 ] 

a = _y t : 

i) 2 2 

2 R 2 W Z 



(2.23) 



If can be expressed as a sum of first and second order polynomials 
in £ , Tj , and t , ( 2 . 22 ) can be integrated directly for all values 
of rcr. r 

The intensity of scattered sound is often expressed with the 
aid of a scattering function S as follows: 




(2.24 ) 



* r 2 2 i 

In Section 6, I replace exp [_- a (l-^)J by the characteristic function 

and remove the restriction that the distribution be Gaussian. 
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where A = illuminated area 



To change (2.22) into this form I need to express A in terms of the 
transducer dimensions and . Since transducers are often described 
by their Fraunhofer "beam width", I do so here. For plane waves, 
defining 0 and X as shown on Figure 2.2, and ignoring time dependence, 





dx' 



dy 1 



(2.25) 



| pj = exp^- k 2 L 2 sin 2 0/4 - k 2 W 2 sin X/ 4 J 

The response is shown on Figure 2.3. I choose to use |p|= e ^ to define 
A0 and AX, as follows 



at X = 0, (kLsin0)/2 = 1 
and at 0=0, (kWsinX)/2 = 1 



Thus sin0 ~ A0 


= 2/ 


(kL) 






sinX ~ AX 


= 2/ 


(kW) 




(2.26) 


The substitution of (2.26) 


into 


( 2 . 22 ) and 


casting into the 


form 


(2.24) yields 










(t + t) 


P (t)) 

CO 


B 2 


A 5 


(2.27) 


D 2 
R 2 





S = k^F ^ 2 e J' J exp j^-a ^ 2 - a^T } 2 + 2 iCC^ + 2 i£r| 

■°° 2 2 “1 

- 4r a (l-! r )J d£ d*rj (2.28) 

A = R^ 2 A0AX/ (cos0^) (2.29) 

Eq. (2.22) is equivalent to Eq. (2.27) and (2.28). We can proceed in several 
ways, direct evaluation at 2.28, by numerical integration, polynomial ex- 
pression of ^ and integration, and special functions for ^ and integration. 
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Fig. 2.2 



Transducer. The transducer is Gaussian shaded. L and W 
indicate the dimensions for e ^ amplitude shading. 
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Fig. 2 




Response of the Gaussian shaded transducer. The far field 
(Fraunhofer) response of the Gaussian shaded transducer, 



g = exp 



jj-x' 2 /L^j 
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Before we cast off, sail into the confused sea, and expand the 
correlation functions, it will pay us to look at the convergence of 



(2.22). The contributions to the integral are small for values of 
i and T) larger than 



6 f > a i 



-% 



T) 



-% 



(2.30) 



and for = R 



0 1 - 0 2 



where 



a^ 2 ~ L/ cos0 






W 



Since the contributions are small, we need not worry about the shape 
of t at distances larger than t and in evaluating the integral. 

Or, the dimensions of the transducer and 0, determine the sensitivity 
of the scattered sound to On the other hand^ accurate fits at small I 
tj are extremely important at large y *a. 



I would like to close this section by giving my thoughts on the 
physical significance of the £, r\ integral for S. This integral 
ought not be confused with the first integral over the illuminated 
area because the I , r| integral relates the phase of the scattered 
signal at any point on the surface relative to a nearby point at 
a displacement I and rj. The contributions to S are small for 
I and tj greater than and t^. Paraphased, the integral is 
in correlation space. For a given roughness, the constancy of the 
phases of the scattered components of signal depends upon the dimen- 
sions of the transducer. 



and 
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Ill- evaluation of the scattering integral 



Experimental measurements of the I, t dependence of wind blown 
water waves have shown that has the form of traveling wave packet, 
Fig. 3.1. The envelope moves at the group velocity and the phases 




Fig. 3.1 * (6, t) 

On the graphs R(t) is our ^ 

of the damped oscillation move at the phase velocity. Near the origin 
(£ > t etc), the dependence on I and t can be approximated as 

♦ 2: *(£ - vt) (3.1) 

for waves traveling in the + x direction, is symetric about 

£ =0, and t = 0 

♦ (£) = >K-6) 

♦ (t) = t(-T) (3.2) 
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Near I = 0, t can be approximated as a polynomial 



2 

♦ ~ 1 * a x - b ( x )[ 

or a*o^T - b|i - 'O'tJ - a^v^T^ (3.3) 

and the expansion has cross terms It. Correspondingly, waves moving 
in the y direction have tjt terms and waves moving in an arbitrary 
direction have both. This doesn't make the analysis more difficult 
because t is a parameter. 

In an earlier paper, I used polynomial fits to the correlation 
function for a random surface (Clay, 1971). The procedure was to 
divide the surface into sub areas and to integrate the scattering 
function for each of the sub areas. I will do the same here except 
that the t dimension is added to the problem. The I, r] and t 
nomenclatures are shown on Figure 3.2. A polynominal expansion for the 
ijkth sub region is 



♦ 



ijk 



. - a ’...I 2 - - e'..^ 2 

ljk ijk ijk ijk 



f ' • -1 T] - g' . t - h 1 T - m' . I t 
j jk ijk ijk ijk 



2 2 



2 2 



• T >’ijk T1T • r ’ijk 5 T - s ijki T 



(3.5) 



It isn't necessary to expand the t dependence as a polynomial in 
some problems. For example, the correlation function can be expanded 
as the product of polynomials in I and functions of t such as cos pT 
and sin pT. If I had the polynomial form programed, I wouldn't bother. 

The coefficients of the variables I and r] can be combined 
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Fig. 3.2 sub areas 

The I , t| map is at constant t. The I, t map is at constant t) . 

V is assumed to be 0 for regions beyond I , ri and t 

max max uidx 
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c = c’ -e 1 t - h 1 t 
ijk ijk 8 ijk ijk 



k,* Air = k * 



. .. + r ' 
ijk 


ijk 


+ m 




ijk 


ijk 


+ s 1 




ijk 


ijk 


. + n' 




ijk 


ijk 



V. .. = c. - a. .A - b ..A 

ijk ijk ijk ijk 



- e . . 



ijk" - £ ijk" 



On designating for the kth t region, we can write 



E 



i J 



. 2 2 2 -iavr 
S. = k % F e 

*-J k 8n 



ijk 
Tlj 



a.r - a V 
5 TJ 



//“ p [• 

6 i-l’ ^J-l 

+ 2iCg + 2ipr] - 4r 2 o 2 (l - '•'ijk)] d ^ dT l 



(3.6) 



(3.7) 



(3.8) 



The integral has the form of the product of two complex 
error integrals. To use the tabulations, we transform the variables. 
Since the integrals on £ and t] have the same form, we use an 
integral on x and let it stand for either. Note that erf(z) 
and erfc(z) are defined by Abramowitz and Stegun(1965) as follows: 



erf (z) 



_ 2 r -t : 

= 'i i e 



dt 
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/ -f' 

e 



dt 



“Z 



w(z) = e erfc( - iz) (3.9) 

We now cast S into the form of (3*9) by defining the constant C 

xjk ° ljk 



u..,, 

ijk 


V. ., >as 
ijk 


follows 


s 


. 2 2 2 
_ k & F 


c 


ijk 


32 


ijk 


C ijk 


= exp £ - 


2 2 

4 T o (1 


U ijk 


J 


f 5i 

f exp 

ti-l 



U 



ijk^l 



xijk 



xx j 



(3.10) 

(3.11) 

dt (3.12) 



where with the aid of (2.23) and (3.7) 
1 



A 2. .. 
x xjk 



2 2 

= at + 4r o a. 

ijk 



2 2 



B . .. = - iCX + 4 r 0 b. 

xxjk ljk 



(3.13) 



Similar expressions for the t) dependence are 



- ■ a r 

tjk 4 ? j 



exp 



" -JL 



Vl L Vjk 



B . r\ 
yxjk 



dr) (3.14) 



2 2 

1 = a + 4r a e . 

— — ' n iJ k 

A yijk 



B . = - ip + 4 r 2 o 2 f . 

yijk ljk 



(3.15) 



Eq. (3,13) and (3.14) have the same form so we drop subscripts and 
cast it into the form of (3.9) by completing the square and changing 
variables as follows: 
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t 



X 



4* s 



A 

s = AB/2 



g(i-l) = .W + s 




A 


(3.16) 


g(i) = _i + s 
A 




TT A s2 

U = Ae 1 e dt 




J 

g(i-D 


(3.17) 


2 

U = A c S erfc [g (i-1)] - erfc [ g(i) ] 


(3.18) 



If g is real, (3 • 18) can be used for the evaluation of U and V. 

This would be expected for the specular direction when (X and P 
are zero. For complex g, we change the form of erfc(g) to w(z) 
in (3.9) . 

Let z = i g (i-1) 

2 

then e S erfc [g (i-1)] = exp [s^ - g^(i-l)] w[ig(i-l)] (3.19) 

Eq. 3.18 with the aid of (3.19) is 

U = A e S / e S ^ w[ ig(i-l) ] - e 8 ^w[ig(i)] 3 

(3.20) 

V has the same form as U. 

The substitution of A . , B . , etc. into (3.16) is simple and 

XI J K XI J K 

the evaluation (3.20) for and V follows directly. I don't 

see much reason to write the expression because there is a chain of 
substitutions through to (3.5) for the t dependence* The analytic 
evaluation of (3.8) is complete. From all of this mess, I would 
not expect the time dependence of the scattered signals to have a 

simple relationship to the time dependence of the surface. 
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IV. COHERENTLY SCATTERED SIGNAL 



For my estimate of the coherently scattered signal, I will use the 
Eckart procedure and calculate the ensemble average of p. This method is 
very general and is applicable to a much wider range of conditions than is 
superfically apparent. By considering the more general problem, it is 
possible to use the dependence of the coherently reflected signal on y 
to determine the probability density function of the surface W. 

I begin by assuming the coherent signal is an ensemble average of 
many transmissions • The sound pressure is given by Helmholtz 

integral and the expansions of R and R 1 , Eq. (2.11 and 2.12). For a 
moderately directive source and £ very small relative to R and R* , I can 
write 



R + R' - R(£ = 0) + R' (^ = 0)+ 2 y^ 



(4.1) 



I use [space] to represent all of the non-random part of x, y, and t 
dependence of the Helmholtz integral, 



and 




r r 



00 



p = [space] ds 

~'-00 



for £ = 0 



(4.2) 



(4.3) 



The latter expression is exact for all configurations and (4.2) is only 
restricted by the approximation, (4.1). The ensemble average of (4.2) is 



p 00 

^p^ = /[space] ds ) 

’-''-oo 



(4.4) 
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The average in (4.4) is given by 







-00 



(4.5) 



and this is the characteristic function of . Hence, we can write the 
following relations 



[< P ) /p o] 






J-J 



<•? 



(4.6) 



W, 




(4.7) 



Since (p^/p^ is a function of positive ycr, Eq. 4.7 should be regarded 

as the sine and cosine integrals from 0 to The normalized coherent 

reflection is the Fourier transformation of the pdf of the surface. 

For a test, I use the experimental data of Mayo, Wright, and Medwin(70). 

p . n 2 2 2 

They measured the surface distribution W^ and U Py> /p I versus 4y a = g, 
where g is often referred to as the roughness parameter. If W^ is Gaussian 
then 



[< P ) /p o] 2 = exp(-g) for Gaussian W^ (4.8) 

however, their data did not fit (4.8). The data agreed quite well with 
(4.8) for g less than 4 (or yo < 1). At larger yo 9 the coherent com- 
ponent is orders of magnitude larger than predicted by (4.8). All aspects 
of the measurements were painstakingly re-examined, including finite 
illuminated area, but no "mistakes" were identified. The non-fit became 
the motivation of this research. For simplicity in a numerical test, I 
approximated their measured distribution function by the linear segments 
shown on Figure 4.1. The substitution and evaluation of (4.6) for the 
piece-wise linear function is routine and omitted here. A comparison of my 
calculation of |{p)/p Q "j ^ and their data are shown on Figure 4.2. The phase 
of the coherently reflected signal probably changes as yo increases. The 
phase is needed in the evaluation of Eq. 4.7. Since the phase is needed 

for the inverse transform, I have not attempted to do it. 
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The procedure we have described, measure the reflected signal at a = 0 
and then the coherently reflected signal for the roughened surface is easy 
to do in the laboratory. Doing this for a rough sea surface would require 
highly accurate measurements and calculations. I suggest that an alter- 
native procedure might be to measure the cross correlation of signals 
at a pair of separated hydrophones. 



\ 




Fig. 4.1. Experimental probability density function of the model sea surface 
a = .45 cm. The solid line is the piece-wise linear 
approximation to the distribution function. Data furnished by 
Mayo, Wright and Medwin. 
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Fig. 4.2. Coherently reflected signal. The data points are from Mayo, 
Wright and Medwin (1970). The solid line is calculated for 
the probability density function shown on Figure 4.1. The data 
are for angle of incidence =45° and a = .45 cm. 
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V. COSINE CORRUGATED SURFACE 



The surface of water waves often looks like a cosine corrugated surface 
for relatively long times and over fairly large areas. If one watches 
carefully^ the phases and amplitudes of the wave change randomly. It is easy 
to make a cosine surface random analytically by letting the phase change 
randomly between each observation (or signal transmission). In addition, 
an ensemble can be formed of waves having different amplitudes. I assume 
the m ^ wave surface is given 

t = t cos K (x - vt - x ) (5.1) 

3 3 m m 



Where v is the wave velocity, is the amplitude, and x^ is the phase. 
The procedure for calculating the scattered sound signal is to substitute 
(5.1) into (2.16) and to expand the result with the aid of the following 
expansion in Bessel functions 



la cos bx T / \ i o 
e = J (a) + 2 

o 



Z <1>nj n 



(a) cos (nbx) 



(5.2) 



The result of integration and manipulation is 

2 r 2 °° 

P = P 0 e ’ P /3y | J o (2r ^m )e ' a /3X + Y (i)T1 J n (2r ^m ) exp 



+ ) (i) J (2y£ ) exp 
Lt n m 

1 



( 2Q - nK) + inK(vt + ) 

4a m 

x 









. lnK(vt + x > 

4a m 

x 



(5.3) 



where 



A ft B e" 1 ^ " k ^ R l + R 2^J 
_ _o 

P o R^ + R 2 



(5.4) 
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O (j >C p- 




(5.5) 



Before I wipe out most of the terms in (5.3) I should discuss them. The 
traveling water wave introduced an infinite series of time-dependent terms 
having frequencies nKv. These terms contribute a modulation of 
the reflected signal. The higher harmonics are more important at larger 
roughness, i.e. If the reflected signal were processed by means of a 

spectrum analyzer, the spectrum would be like that shown on Figure 5.1. 

The powers in the components depend on the bistatic geometry. 



• O 9 



1 



(co-a^ir) Oo-Ktrj (uu+Ktr) + 



• O * 



frequency 

Fig. 5.1. Spectrum of the signal reflected at a traveling water wave. 

to is the frequency of the sound signal and Kv is the frequency 
of the water wave. 
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The simplest ensemble average is obtained by letting all phases of the 
traveling wave be equally likely. An average over Kx^ from 0 to 2 it eliminates 
all of the terms in the summation in Eq. (5.3). In the specular direction 
(5.3) then reduces to 



(p) /p n = J (2rt ) 

\ / o o m 



(5.6) 



The distribution function of a cosine wave is 



t = t cos K x 
b b m 



W 

c 



-1 






i < 5 



m 



W = 0 
c 



i > i 



m 



(5.7) 



The circle can be completed by averaging (5.6) over a random set of 
amplitudes. For example, the envelope of a Gaussian random function has 
a Rayleigh distribution function: 

W R = o~ 2 5 exp [- £ 2 /(2 o 2 )] (5.8) 

f J (2y{ ) W R d£ = exp j^- 2y 2 a 2 J (5.9) 

J o ° 

This would be expected on the basis of the central limit theorem. I suggest 

that fair approximations to non-Gaussian functions can be made by averaging 

(5.6) for several values of t . 

3 m 
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VI. TOTAL SIGNAL SCATTERED AT A NON-GAUSSIAN SURFACE 



This is my last section and lest the reader has forgotten, my purpose 
is to do the inverse problem. The basic formulation of the scattering 
problem is given in Sections 2 and 3. There, the emphasis was on going 
from a known surface to calculated scattered signal. This analysis 
started with a bivarate Gaussian surface and its correlation function. 

I could have used the characteristic function and removed the Gaussian 
restriction as follows: 



c 2E ^ 21 r(C-0> «. 

The only change is to replace exp£ -4y 2 a 2 (l-^)3 by in Eq. (2.22). 
Following my procedure in Section 4, I evaluate p^ and express Eq. (2.22) 
as follows: 

/ A / 2 , sh -1 -iayr 

\PP'y /P 0 = (a^ a ^) 2 it e 






exp 




c d ? d 
2 £ T) 



( 6 . 2 ) 



D = exp [- a^l 2 - a n 2 ] 



(6.3) 



S' T, 

Since a, and a are functions of k, , 9^, and 9^, they can also be 
5 T) 1 2 3 J 

expressed as a function of CL and (3. I make a Fourier transform operation 
on assuming that D and D ^ are slowly varying functions of i and *q. 
Integration of CL and f3 yields 5(£ - £') and 5 (t] - tj ' ) type functions and 
the expression for is the following 



— L_ rr P( a a iLJT r . , , , . , ,-i 

2k 2JJJ ° 2 e exp L' l(cj T + 2o ^ + 2p-n)J dadpd 



cj * = C, 



(6.4) 
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The G! and p dependences of a^ and a^ are greatly simplified by letting 
9^ = ± ^ 2 > ^3 = and folding and constant for a set of measurements. 
Obviously, one * s ability to do the inverse transformation depends upon 
having accurate values of ^pp*^ . I assume the user would do the trans- 
forms numerically. Although the integral has infinite limits, actual 
measurements of ^pp^ give both upper and lower limits. There is another 
consideration. Since D operates as a spatial high pass filter in (6.2), 

I doubt if would have much accuracy at I and T) much larger than and 

v 
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